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Homogenization protocols model the quantum mechanical evolution of a system to a fixed state 
independently from its initial configuration by repeatedly coupling it with a collection of identical 
ancillas. Here we analyze these protocols within the formalism of "relaxing" channels providing an 
easy to check sufficient condition for homogenization. In this context we describe mediated homog- 
enization schemes where a network of connected qudits relaxes to a fixed state by only partially 
interacting with a bath. We also study configurations which allow us to introduce entanglement 
among the elements of the network. Finally we analyze the effect of having competitive configura- 
tions with two different baths and we prove the convergence to dynamical equilibrium for Heisenberg 
chains. 

PACS numbers: 



I. INTRODUCTION 

Homogenization protocols have been extensively stud- 
ied in recent years as a powerful model for the equilibra- 
tion of a quantum mechanical system interacting with a 
large bath [J H, H, 0, H, El ■ In these schemes one considers 
a collision-like coupling of the system with a collection 
of ancillary systems that have been prepared in iden- 
tical states. This corresponds to a Markovian approx- 
imation in a discrete dynamical evolution. Compared 
to typical quantum Markov equations the advantage of 
this model is that it allows one to concentrate on the 
effective unitaries and completely positive maps rather 
than the underlying Hamiltonians and Lindblad genera- 
tors Q- By "homogenization" one means that the sys- 
tem state converges to a state that is the same as the 
ancilla states. This was demonstrated for a class of qudit 
systems in 1 , 0, H[ . The bath-system entanglement was 
studied in J], and a continuous-time model (quantum 
master equation) was derived from the discrete model 
in Furthermore, the emergence of irreversibility was 
investigated in Q. 

An important aspect of quantum homogenization is 
that is as a stable method of driving a system into some 
fixed state, independent of its initial state. In this con- 
text, we will also refer to the bath as a "controller" sys- 
tem and its state as a "controller state". Hence apart 
from its fundamental role of studying quantum con- 
vergence, homogenization has possible applications for 
quantum cloning [2| , for the hiding of quantum informa- 
tion 0, [f| and spin chain quantum communication Q . 

The prototypical homogenization scenario [J, 0] is de- 
scribed in Fig.[TJ It is composed of two parts: a system A 
with an always on Hamiltonian Ha, and a large ensem- 
ble of identical controller systems B%, B2, ■ ■ ■ , B n . The 
latter are prepared in the same initial state u>b and are 
assumed to have no independent free evolution. The sys- 
tem A is coupled in sequential order with each one of 
the Bs through a series of identical stepwise interactions 
described by the Hamiltonian Hj. In this setting the 
evolution of the system A is described by the successive 
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FIG. 1: Setup of a standard homogenization protocol: the 
controlled system A interact with a collection of controller 
systems B which have been initialized into the same input 
state u>b- Homogenization takes place when the final state of 
A is driven into the same state u A of the controllers in the 
limit of infinitely many couplings with the Bs. 



application of the completely positive (CP) map 

£(j>A) = T±B[U(pA®w B )Ui] , (1) 

with U = cxp[— i(H A + Hj)t] and t > being the time 
interval associated with a single A-B coupling. After the 
interaction with n controllers the state of A becomes 



=£o£o-_o£( PA ) = £ n {p A ) 
n times 



(2) 



We are interested in the behavior of the sequence ([2]) 
for large n: in the case where the system A and the con- 
trollers Bk are identical, Ha = and Hj is a swap Hamil- 
tonian it was shown 

0, S 3 that the state of A asymp- 
totically converges to the state wg of the controllers, in- 
dependently from the initial state pa- 

In the above, homogenization also gives rise to ther- 
malization - if the bath is initialized in Gibbs states, 
then the system converges to a Gibbs state. However 
it is an open question if this still holds in a situation 
where system and bath particles have different dimen- 
sionality [||. Moreover, in [J, H, i, i, i, i] the bath is 
modeled to interact with the whole system whereas an 
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FIG. 2: Generalized homogenization protocol: in this case 
the controlled system A is a composite one (e.g. a network of 
coupled spins) . Only a proper subset of the network interacts 
directly with the controllers B (as in the case of Fig. [TJ the 
Bk are assumed to be prepared in the same initial states.) 



interaction with a subsystem (such as the surface) seems 
more plausible. 

A first generalization towards this direction was ob- 
served by the Authors of the present paper when study- 
ing the propagation of quantum information along spin 
chain communication channels In that case A repre- 
sented a collection of N coupled qubits, while the Hamil- 
tonians Hj implemented a sequence of strong instanta- 
neous swaps among the last element of the chain and 
a collection of controller qubits (the Bs). By assuming 
the Bs to be prepared into the spin down state |0)s we 
showed that in the limit of large n, any initial state of 
A will be coherently transferred into the Bs while the 
chain will be mapped into the all spin down configura- 
tion 1 00- • • 0) (the only requirement being a non trivial 
connection among the qubits of A) . Concerning Ref . M it 
is worth stressing that in contrast to Refs. [H, H d, 0, H Q 
A and the controllers are quite distinct objects (namely A 
is a network of coupled qubits while each of the B\ ■ ■ ■ B n 
is just a single qubit). Moreover only a proper subset of 
A (specifically the N-th. element of the network) interacts 
directly with B: the remaining qubits are only affected 
by the controllers through the free Hamiltonian Ha of 
the network (see Fig. [2]) . The possibility of preparing the 
ancilla state |0) b into each one of the spins of the network 
is therefore a remarkable feature of the system which re- 
quires some further investigation. We call it mediated 
homogenization process. 

In this paper we tackle this issue by analyzing spin 
networks which show similar properties. In particular 
in Sec. IIIII we present a first example of a mediated ho- 
mogenization process which allows one to transfer on the 
network element any input state lob of the controllers. 
Before doing so however we introduce a general conver- 
gence criterion for relaxing channels in Sec. [TT] that will 
turn out to be extremely useful in our discussion. Our 
generalized setup gives rise to much richer quantum con- 
vergence effect - in Sect. IIVI we give numerical evidence 
for equilibrium states which are entangled. Finally in 
Sect. |V] we study the effects of having competitive baths 
at different temperature. 



II. MIXING CRITERIA IN THE 
HOMOGENIZATION SETUP 

Our starting point is the observation that what makes 
the CP map in Eq. @ converge into a specific state in- 
dependently from the initial input pa is a well known 
property called relaxing [lOj (also referred to as "mix- 
ing" [ll|, [l2j or "absorbing" [lj] ) . In this language the 
convergence point 

Pa = li m P { a > ( 3 ) 

n — >oo 

is the only fixed point of £, i.e. the only solution of the 
equation 

£(j>a) = pa , (4) 

(we refer the reader to Refs. [ill, EH for an detailed in- 
troduction to relaxing channels). Therefore an homoge- 
nization procedure is associated with a relaxing channel 
whose fixed point uja coincides with the state of the con- 
troller. 

We now prove a very simple but important result. Sup- 
pose that given U we know that the map ([I]) is relaxing 
for a specific choice of the controlled state lob- Consider 
now the map £ which is obtained from Eq. (JTJ) by replac- 
ing ujb with a state ujb which is a (non trivial) convex 
combination of pb, i.e. ljb = P + (1 — p) w' B , with 
p S]0, 1] and lo' b being a generic density matrix. In this 
case the map £ can be expressed as a convex combination 
of £, i.e. 

£ =p£ + (1-p)£', (5) 

with the map £' as in Eq. {TJ) with lob replaced by lo' b . 
We can then use a theorem by Haag [lj] which shows that 
a convex combination of CP maps containing at least one 
relaxing channel is relaxing, to conclude that £ is relaxing 
— for completeness, we provide an alternative (and much 
simpler) proof of this important theorem. It is based 
on the fact relaxing map is equivalent to the asymptotic 
deformation property Hence for all p ^ p' there is 
a k such that 

||£ fc (p)-£V)lli<llp-p'lli, (6) 

where ||0||i = Tr[V 0t0] indicates the trace norm of the 
operator 0. We write £ k — p k £ k + (1 — p k )S' , where 
S' is the CP map that contains all other terms of the 
expansion of £ k . By the non-expansiveness [15( of S' and 
the triangle inequality we obtain 

||f k (p)-fV)IU<IIP-p'lli, (7) 

whence £ is an asymptotic deformation. In the context of 
Fig. [2] this implies that the system still converges when 
substituting the controller state lob with ujb- For in- 
stance, if there exists any state lob for which the system 
converges, then it will converge to the fully mixed state 
if the Bk are initialized in the fully mixed state. 
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A natural question is then to determine the fixed point 
of £. Specifically one may ask how the final state of the 
system A depends upon the controller state oib- For 
instance: if homogenization takes place for oib, does it 
hold also for cDg? Or, how does the entropy of the fixed 
point depend on the entropy S(o>b) of the controllers? 

Before passing to apply the Haag criterion to the me- 
diated homogenization scheme, it is worth presenting yet 
another interesting generalization of this simple but im- 
portant theorem. Consider in fact the situation in which 
the states of the controllers Bi , • • • , B n have not being 
properly initialized. In particular we are interested in 
studying what happens if instead of being prepared in the 
"good" initial state oib, the £th controller is described by 
the following imperfect state 



01$ = pe oj b + (1 -pi) 



Qb 



(8) 

(£) 

where for £ = 1, • • • ,n, pi > are probabilities and g B 
are density matrices. According to the analysis of Sec. U 
this yields a sequence £\ ■ ■ ■ , £„ of CP maps which have 
the property that each of them is a convex combination 
of a fixed relaxing map £, i.e. 

El = Pt S + (1 - Pt )St , (9) 

where £ and Sg are respectively the channels ([1]) asso- 

(£) 

ciated with ojb and g B respectively. Clearly without 
putting any restriction on the values of pi nothing can 
be said about the convergence property of the proto- 
col. Therefore we consider the case in which the "error" 
(1 — pi) is bounded, by imposing the constraint 



pi ^ p > 



(10) 



This hypothesis does not yet guarantee that A will be 
driven toward oi A - However we can at least verify that 
the process is still able to "forget" about the initial state 
of the controlled system as in the relaxing case (this is 
a typical feature of any homogenization protocol). The 
evolution of the controlled system A is in fact now de- 
scribed by the following sequence of concatenated maps, 



Mr 



En E n -1 O ■ ■ ■ o£ 1 



(11) 



For arbitrary input density matrices p' A , p" A of the con- 
trolled system define 



fn = \\M n {p' A )-M n {p' A )\\ 



(12) 



Now since /„ is non-negative and non-increasing [15| it 
certainly admits a limit lim TWOO /„ = /». To show that 
the protocol forces the controlled system to forget about 
its initial conditions we need only to verify that this quan- 
tity is null for all p' A and p" A . Assume then by contradic- 
tion that /* > for some choice of these input states. Let 
p* A be the fixed point of the unperturbed map £. Then 
there is a value of k such that \\£ k (pA) — p A \\ i < /*/4 for 
all pa OH- Let then S = 3 ^J* k ^ > 0. There is a n such 
that /« — /*< 5. We write 

M k+ n = MoM n (13) 



where the superoperator M = £k+ n ° 
decomposed as 



■ o £i +n can be 



M = P k £ k + (1 - P k )T (14) 

with T being CP and P k = p k + n ■ ■ ■ Pi+ n > p k by as- 
sumption. Hence 



f k +n = \\M(Mn(p A )) ~ MiMnU 



A)\\l 



< p k \\e k (.M n (PA))-e k (M n (// A ))\\i 

+(1 - P k )\\T(M n ( PA )) - T(M n (p A 

< P k f*/2+{l-P k )(S + f*) 

P k f* 



< ft/./2+(l-flfc) 
= f* ~ Pkf*/6 < U 



3(1 - P k ) 



f* 



Since /„ is non-increasing this is a contradiction, and 
/* = 0. We have shown that the whole state space is con- 
tracted to a single point. In general, this point is still 
evolving under the action of £ ni but contains no infor- 
mation about the initial state. The map M n is relaxing 
if and only if there exists an asymptotic fixed point g A , 
i.e. a state with lim^oo A4 n (g A ) = g* A . 



III. MEDIATED HOMOGENIZATION IN SPIN 
NETWORKS 

An interesting example of mediated homogenization is 
obtained by assuming A to be a network of N coupled 
qudits A\, - ■ ■ , An mutually interacting through a sum 
of local term of the form 



Ha — ^ J kk 'SA k 



(15) 



fc.fc' 



where J kk > are coupling constants and where S Ah A k , — 
(SA k A k , )' are unitary operators which swap the k-th qu- 
dits of A with the fc'-th [16[ . Regarding the coupling with 
the controller we consider the case in which only the An 
interact with the Bs (also represented by d-dimensional 
systems) through a swap Hamiltonian similar to (|15p . i.e. 



Hi — Sba* 



(16) 



Under these conditions we can show that, for all choice 
of the controller states o>b and for almost all choices of 
the interaction time t the map (JT|) is relaxing with fixed 
point 



P*A = {"A) 



(17) 



given that the graph associated with the coupling J kk i 
satisfies certain constraints. This corresponds to the case 
in which, in the limit of large n, the controller state oib 
is "copied" in all the TV controlled qudits. We call this 
process a mediated homogenization of A. It fulfills all four 
homogenization criteria mentioned in Ref. 6]: Firstly, 
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the coupling between system and bath is independent of 
the bath state. Secondly, the equilibrium state is not only 
a fixed point of the CP map |T]) but also of the unitary 
evolution?/ = cxp[— i(Ha + Hj)t] alone. Thirdly, the 
system converges to the fixed point for all initial states. 
Finally the change of the bath due to the evolution can 
be made arbitrarily small by choosing a short interaction 
time t. An immediate consequence of the above result is 
the fact that the von Neumann entropy of A converges 
to N times the von Neumann entropy of the controller 
state Sb = — Tr[/9B log 2 Pb\- This is a distinctive trait of 
the mediated homogenization processes and it is similar 
to what happens when we put a system of interest in 
thermal contact with an reservoir. It should be pointed 
out though that in our case the convergence state is in 
general far away from any thermal state exp[— [3Ha]/Z. 
The thermalization feature observed in [l[ thus seems to 
be specific to the case where A is a single qudit. 

To prove the above result we first focus on the case in 
which lub is a pure vector \4>)b- Define then the joint 
observable 

M AB = M B + M A , (18) 
where Ma — J2k=i ^A k and 

M B = -\<p) B {<t>\ 

M Ak ee -\4>) Ak {4>\ . (19) 

The operator Mab commutes with the total Hamilto- 
nian H = Ha + Hbc and hence with the operator 
U = cxp[— iHt]: we thus say that free evolution of the 
network preserves the "excitations" associated with the 
projectors \<f>){4>\- 

Moreover the state \<j)) b is the (non degenerate) eigen- 
vector associated with the minimum eigenvalue of Mb- 
Under these conditions we can invoke the Lemma 3 
of Ref. [ll| which states that the map |T|) is relaxing 
with fixed point \<ft)® N = \<j>) Al <8> ■ • • ® \4>)a n if the 
state \(j>)® N (g) \4>) b is the unique eigenvector of U hav- 
ing the form \E)a ® \4>)b- This last condition can be 
verified by focusing on the global Hamiltonian H: if 
indeed \</>)® N <£> \4>)b is the unique eigenvector of the 
global Hamiltonian H with the factorization property 
\E)a® |</>)b then the same property will holds for U for 
almost all the values of t. We have shown elsewhere [TtJ 
that for "excitation" preserving Hamiltonians the above 
factorization condition depends only on the geometry of 
the associated graph. For example, an open chain with C 
being an end qudit has the required property. We can ap- 
ply this result to the Hamiltonian (TT5|) ( it is excitations 
preserving in the sense that it does preserve the number 
of qudits of A which are in the state | </>)). Therefore for 
any given network configuration satisfying the topologi- 
cal constraint of Ref. [17| we can conclude that, for all 
initial pure state |</>) of the controller, the above iterative 
procedure will drive A to a unique fixed point. Before de- 
termining such fixed point, let first observe that the same 
convergence will hold also when assuming the initial state 



of the controllers to be a general mixed state u>b- This is 
a trivial consequence of the pure case scenario which can 
be obtained by expanding any such mixture into a convo- 
lution of pure states u>b = J2j Qj^j) B{4>j\ and applying 
the Haag criterion. 

Since we have now proved that for all choice of the 
controller state lob the channel £ is relaxing, to verify 
Eq. (fl~7|) it is sufficient to show that (lja)® N satisfies 
Eq. ([5]). The latter can be easily verified by noticing 
that each summand of the Hamiltonian (|15p commutes 
with all tensor product operators of the form <d® N , and 
therefore 

[H A ,Q 9N ]=0. (20) 

Consequently for p* A as in Eq. (fT7)l and Hbc as in 
Eq. (|16|) we can write 

[H A + H BC , P*a®Pb] = 0=> [U, p* A ®p B ]=0, (21) 

which is sufficient to show that p* A satisfies the invariancc 
condition (g]). 

IV. BUILDING ENTANGLEMENT IN THE 
NETWORK 

In the previous section we found a model where inde- 
pendently from the initial state of A, the final state of the 
network is the separable state w® N . Each of the N qu- 
dits of the network has been driven into the initial state 
of the controllers. In this section we show that, keeping 
Hj as in Eq. (fTB)) , there are also Hamiltonians H A which 
are capable of building entanglement among the qudits of 
the network. Although this is no longer a homogeniza- 
tion protocol (the controller state is not transferred to 
the controlled system) it could have useful applications 
as a method of state preparation. 

Consider for the sake of simplicity d = 2. In this case 
the swap interaction of Eq. (|15p corresponds (up to a 
constant) to a Heisenberg coupling. A natural generaliza- 
tion of it is then provided by the anisotropic Heisenberg 
Hamiltonian 

ha = y, Jj y- {"P'P + *W + A ^) ' (22) 

k,k' 

where cr^ ,y ' z ^ represents the Pauli matrix of the A:-th 
qubit of A and where A — 1 is the anisotropy param- 
eter (the isotropic coupling is obtained for A = 1). For 
this coupling we can use the same argument given in 
previous section to characterize the relaxing properties 
of the associated map ([1]) (in particular the factorization 
property of its eigenvectors depends only on the geom- 
etry of the associated graph |17j). In this case however 
the isotropy is lost and the Hamiltonian has a preferred 
spatial direction associated with the z axis which makes 
|0)b and \1)b special with respect to the other controller 
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pure states. Indeed we can still show that mediated ho- 
mogenization takes place for input states uj b which are 
diagonal in the computational basis, i.e. 



WB =p|0>B<0| + (l-p)|l>B<l| 



(23) 



This follows by the fact that for such a choice Eq. (j2Tj) 
holds independently from the value of A. On the contrary 
for more general controller states mediated homogeniza- 
tion is lost. As an example, consider 



uj B =p\0) B (0\ + (l-p)\-) B (-\ 



(24) 



where \—)g = (|0) b — |1) b) /V% and p > 0. In this case 
Haag's theorem can still be used to ensure relaxing of 
the map (|TJ) even though computing the fixed point is 
not simple. For such choice however we have numerically 
verified that the mediated homogenization does not take 
place in general. We evaluated the asymptotic limit of 
the von Neumann entropy of p* A , verifying that it is no 
longer a multiple of the von Neumann entropy of the 
bath state. In Fig. [3] we show an example for a XX 
chain (A = 0). 

Of particular interest is the case p = 0. In this limit 
lub = \ — )b(~ I and the relaxing property cannot be es- 
tablished from the theorem by Haag (simply ps is not a 
convex combination of |0)s(0|). Nevertheless we can use 
numerical analysis to show that the map |T]) is still relax- 
ing [10| . We found that the convergence point is highly 
mixed. Since this is so much different from the isotropic 
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FIG. 3: Ratio R of the "output" entropy Sa of the conver- 
gence point and the "input" entropy Sb of the bath state for 
a XX chain of length N = 4. The controller state is as in 
Eq. (|24p . For all values of p the dynamics is relaxing — for 
p > this is a trivial consequence of Haag theorem, for p = 
instead it can be directly proved by numerical means [l8| . 
When p — > 1 the state becomes diagonal and the ratio con- 
verges to N. For p — > the ratio diverges as the bath state 
becomes pure but the convergence point remains mixed. This 
should be compared with the behavior of a mediated homog- 
enization process (e.g. the swap coupling of Sec. IIII|) where 
R is always equal to N for all p. The parameters for the 
numerics are J kk ' = <5fe,fc+l an d t = 0.5. 

Heisenberg model with fixed point |— )® N it seemed nat- 
ural to compute the relaxing property and convergence 
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FIG. 4: The entropy Sa of the convergence point and the 
concurrence between the first two qubits for an anisotropic 
Heisenberg chain of length N = 4 as a function of the 
anisotropy parameter A. The bath state is given by (|0)s — 
|l)s)/v / 2. The parameters for the numerics are J k y = 8 k ^ k+1 
and t = 0.5. 



point of the anisotropic model for p = as a function 
of A to see the transition for a XX chain (with highly 
mixed convergence point) to the Heisenberg chain (pure 
convergence point). In particular we wanted to check if 
there are also entangled fixed points. For this purpose we 
computed the concurrence between the first and second 
qubit of the chain (say) for intermediate A (see Fig. [4|. 
Again, for the given parameters, all examples were re- 
laxing. We found that the convergence point is indeed 
entangled for some values of A. Contrary to the results 
in the last section, the numerical examples of conver- 
gence points observed here depend on the parameters of 
the model. An important open problem is to determine 
if there exist Ha and lob that have a fixed point with 
interesting applications (e.g. a cluster state). 



V. DYNAMICAL EQUILIBRIUM 

The many-body structure of A presented in Fig. [2] al- 
lows us to consider more complicated procedures. For in- 
stance we can analyze competitive configurations where 
the dynamics of the network A is driven by the simulta- 
neous coupling with two independent sets of controllers 
(the Bi, ■ ■ ■ ,B n and the C\, ■ ■ ■ , C n of Fig. [5]). We can 
then model the "transport" of excitations through the 
network by assuming the two sets to be directly coupled 
with distinct network elements (say Am for B and A-y 
for C) and assuming different "temperature" for the two 
species of controllers (say ujb — p|0)s(0| + (1 — p)|l)s(l| 
for the Bs and v c = q\0) c (0\ + (l-g)|lic(l| for the Cs). 
A similar situation is considered in [n| [2(| Hl| where 
in the case of a linear chain coupled through Heisenberg 
and XX interactions the relaxing property was observed 
numerically (22j. Here, the convergence can be derived 
analytically for arbitrary chain length as a consequence 
of the Haag theorem. To verify this it is convenient to 
treat B and C as a unique controller composed by el- 
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FIG. 5: Setup of the dynamical equilibrium: here the sys- 
tem of Fig. [2] is coupled to two competing baths at different 
temperatures. 
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FIG. 6: Probability of finding the fcth qubit in the state |0)(0| 
for the steady state of a Heisenberg spin chain of length TV = 
5. The plot includes the systems B (site 1) and C (site 7). The 
parameters p — 0.9 and q = 0.4. We give two examples with 
different choice of the interaction time t. 

ements B\C\, B2C2, • • • , B n C n . From the above defini- 
tions it then follows that such composite controllers are 
initialized in the state 

u B ®vc= pq \0)b(0\ ® |0)c(0| + (1 - pq) Qbc , (25) 

with qbc being a density matrix. Therefore according to 
the Haag theorem the convergence can be verified by fo- 



cusing only on the case in which B and C are initialized 
in |0}s(0| (g> |0)c(0|. With this choice however the iter- 
ative procedure is equivalent to the "cooling" protocol 
discussed in Ref. [17| and the convergence is automat- 
ically verified. Deriving the exact steady state in this 
case is however quite complicated so we restrict to nu- 
merical analysis. Again its form depends strongly on the 
parameters, as shown in Fig. [5J 



VI. CONCLUSION 

We have generalized the homogenization protocols to 
a scenario where the system is no longer a single qu- 
dit. We found that mediated homogenization still takes 
place on the lattice when the interaction is taken to be 
isotropic. Anisotropic interactions on the other hand do 
not in general show homogenization. Our numerical re- 
sults are quite suggestive in this direction but are cer- 
tainly not conclusive. This suggests many further stud- 
ies: what is the structure of the fixed points of these 
systems? How are their entropies related to the bath 
entropy? What happens when the system is close to a 
critical point? Can we use this convergence as a way of 
preparing useful states such as cluster states on optical 
lattices? Finally we looked at transport along chains in- 
terconnecting baths at different temperature, where the 
Haag criterion allowed us to prove the convergence to a 
dynamical equilibrium. We found that the temperature 
profile is strongly depending on the parameters of the 
system, such as the interaction time, and not even mono- 
tonic for some times. While at the moment these results 
are numerically only, it may be possible to obtain an an- 
alytic expression for the fixed point in a weak coupling 
limit by deriving a closed equation for the proper ansatz 
(cf. [2y|)- This will be subject of future investigations. 
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